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A  Perturbation  Theory  for  Light  Diffraction  from 
a  Bigrating  with  Multiple  Surface-Polariton  Excitation 

N.  E.  Glass 
Department  of  Physics 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Abstract 

A  perturbation  theory  has  recently  been  developed  for 
treating  the  diffration  of  light,  with  surface-polariton  resonant 
excitation,  from  a  bigrating  surface.   The  theory  is  an 
approximation,  to  first  order  in  the  grating  height,  of  an  earlier 
theory  based  on  the  Rayleigh  hypothesis.  The  perturbation  theory 
holds  for  arbitrary  polarization  and  for  arbitrary  plane  of 
incidence  with  respect  to  the  grating  directions.   It  is,  however, 
limited  to  treating  only  the  possibility  of  two  simultaneous 
resonant  evanescent  waves,  whereas  four  are  possible  at  normal 
incidence  on  a  square  grating,  and  three  are  possible  on  a 
rectangular  grating.   The  present  work  generalizes  the  earlier 
work  to  allow  for  a  four-fold  (or  three-fold)  resonance.   This 
extension  also  allows  one  now  to  determine  the  complex  dispersion 
relation  for  surface-polaritons,  for  wavevectors  at  the 
intersection  of  two  Brillouin  zone  boundary  lines  (not  possible  in 
the  previous  formulation) . 


I.  Introduction 

The  study  of  electromagnetic  wave  scattering  from  corrugated 
metal  surfaces,  with  resonant  excitation  of  surface-plasmon- 
polaritons,  has  many  important  technological  applications.   One  of 
these  is  the  creation  of  electromagnetic  radiation  from  an 

electron  beam  passing  above  a  grating  surface  --  the  Smith-Purcell 

(1)  •    . 

effect v   .   In  many  applications,  the  enormous  enhancements  in  the 

absorptance  at  the  surface  and  in  the  field  intensity  in  the 
surface  region,  due  to  the  grating  coupling  of  an  incident  wave 
into  a  surface-polariton  (or  surface-electromagnetic-wave,  SEW) , 
can  be  further  enhanced  by  use  of  a  bigrating  (a  surface  periodic 
in  two  different  directions)  in  place  of  a  classical  grating 

(periodic  in  one  direction) .   It  was  pointed  out  in  some  of  the 

( 2-3 ) 
earlier  theoretical  studies  on  bigratingsv    '  that  the  resonant 

absorption  of  unpolarized  light  (e.g.,  of  sun-light  in  solar 

absorbers)  will  be  greater  on  a  bigrating  than  on  a  classical 

grating.   When  the  plane  of  incidence  is  normal  to  the  groove 

direction,  only  the  p-polarized  component  of  the  incident  wave  can 

excite  the  surface-polariton  (which  is  itself  a  p-polarized 

surface-wave) ;  the  s-component  does  not  couple  to  the  surface-wave 

at  all.   By  contrast,  the  bigrating  can  couple  both  the  p-  and  the 

s-  components  of  the  incident  wave,  regardless  of  the  orientation 

of  the  plane  of  incidence.   Moreover,  as  demonstrated 

( A) 
experimentally  by  Inagaki  et  al.    ,  the  bigrating  can  couple  the 

incident  wave  simultaneously  into  two  surface-polaritons,  leading 


to  an  extra  absorptance,  for  incidence  orientations  in  which  the 
classical  grating  would  only  excite  a  single  polariton.   These 
same  processes  indicate  mechanisms  by  which  the  Smith-Purcell 
radiation  can  also  be  enhanced  by  using  a  bigrating  in  place  of  a 
classical  grating. 

The  first  theory  of  the  Smith-Purcell  effect  that  was  both 

mathematically  exact  and  physically  realistic  was  the  work  of  P.M. 

(5)  . 

van  den  Bergv  '.   This  theory  used  a  Green's  function  approach, 

(  6 ) 

proposed  by  di  Francia^   ,  in  which  the  component  evanescent 
waves,  which  make  up  the  Fourier  expansion  of  the  field  of  the 
electron  current,  are  seen  as  exciting  the  various  diffracted 
orders  of  the  classical  grating. 

One  shortcoming  of  the  van  den  Berg  theory  is  that  it 
considers  the  metal  on  which  the  grating  is  ruled  to  be  a  perfect 
conductor.   This  approximation  precludes  the  possibility  of 

accounting  for  the  effects  of  the  surface-polaritons  and  their 

(7) 
resonant  excitation.  Recently,  Chuang  and  Kongv  '  developed  a 

theory  which  considers  a  sinusoidal  grating  on  silver,  described 

by  a  complex  dielectric  function,  and  uses  an  extended  boundary 

condition  method  to  calculate  the  radiated  energy. 

Similar  techniques  are  in  principle  applicable  to  the 

treatment  of  Smith-Purcell  radiation  from  a  bigrating,  but  the 

numerical  difficulties  can  become  prohibitive:  the  computing  time 

necessary  to  calculate  the  diffracted  field,  arising  from  each 

Fourier  component  in  the  incident  field,  will  be  much  greater  for 

the  two-dimensional  case,  the  bigrating,  than  for  the 


one-dimensional  case,  the  classical  grating.   Thus  perturbative 
technigues  are  desired. 

A  non-perturbative  theory  of  light  diffraction  from  a 

bigrating,  with  surface-polariton  resonances,  was  developed  by 

( 8 ) 
Glass,  Maradudin,  and  Celliv   .   This  theory  uses  the  vectorial 

form  of  Kirchhoff's  integral  to  eliminate  the  field  in  the  metal 

and  thereby  to  cut  in  half  the  size  of  the  matrix  eguation  to  be 

solved.   There  are  nevertheless  very  large  matrix  eguations  to  be 

solved  with  this  method:  from  98x98   up  to  243x243  in  Ref.  8. 

(9) 
Thus,  a  perturbation  theory  was  recently  developed  by  Glass v   , 

which  applies  the  coupled-mode  theory  of  Glass,  Weber,  and 

Mills^   '  to  the  exact  theory  of  Glass,  Maradudin,  and  Celli^   . 

This  perturbation  theory  reduces  the  computations  to  a  4x4  matrix 

eguation,  by  retaining,  to  first-order,  the  two  amplitudes  (p-and 

s-  polarized  components)  of  the  specularly  diffracted  beam  and  two 

evanescent-wave  amplitudes  (p-components  only)  corresponding  to 

the  possible  simultaneous  resonant  excitation  of  two  surface- 

polaritons.   These  two  surface-polaritons  may  be  propagating  in 

non-collinear  directions  (unlike  the  case  for  the  classical 

grating,  in  which  two  polaritons  can  be  excited  only  for  normal 

incidence,  and  then  are  propagating  in  exactly  opposing 

directions).   This  theory  as  implemented  in  Ref. 9  has  one 

shortcoming:  it  cannot  treat  the  case  of  more  than  two 

surface-polaritons  excited  simultaneously.   For  example,  it  cannot 

treat  the  cases  of  normal  incidence  on  a  sguare-grating  when  there 

are  four  surface-polaritons  excited  simultaneously  or  of  obligue 


incidence  on  a  rectangular  grating  when  there  are  three 
simultaneous  excitations.   This  shortcoming  could  be  an  important 
one  if  one  wishes  to  use  these  perturbative  approaches  in 
calculating  the  Smith-Purcell  radiation.   Fortunately,  it  is  not  a 
difficult  problem  to  correct.   A  straight-forward  extension  of  the 
theory  of  Ref.  9  will  produce  a  6x6  matrix  equation  for  four 
resonant-wave  (p-)  amplitudes  and  the  two  specular  beam 
amplitudes.   It  is  the  purpose  of  this  report  to  present  the 
results  of  the  extension  of  the  perturbation  theory  of  Ref.  9  to 
treat  the  cases  in  which  there  are  up  to  four  surface  excitations. 

II.  Theory 

The  physical  geometry  is  as  follows.   The  nominal  flat 
surface  is  in  the  plane   x  =0.   The  surface  profile  is  described 
by   x  =C(Xn)   [where   x(|=x  x  +x  x    ].   The  region   x  >C(x..)   is 
vacuum;   x_<C(x..)   is  dielectric,  described  by  the  complex 
dielectric  function   e (w) =e  (w) +ie  (w) .   The  surface  profile 
function  is  periodic  in  two  directons.   Light  of  frequency  w   and 
wavevector   k   is  incident  from  the  vacuum  (with  angle  of 
incidence   9).   The  projection  of   k   onto  the  surface  is   k|(  , 

A 

which  makes  an  angle   4>   with   x  . 

The  electric  field  in  the  vacuum  is 

f!(x)  =  ^(g),]?!,)  exptik^-Xn  -  iaQ  (co  ,  k„  )  x3  ]   + 

+  |  fS(w,^)  exp[i^-x„  +  iao(u,K£)x3]   .        (1) 
G 

In  the  above  equation  we  have: 


&  «  8..  +  3  ,  (2) 


G 
where   G   is  a  reciprocal  lattice  vector  defined  by  the  bigrating 

periodicity,  which  for  a  simple  square-lattice  grating  is 

G  =  — -  (m  x   +  m  x  )     for  m .=0,±1,±2 , . . .  (3) 

("a"  is  the  period) ;  and  where 

k|(    =   K     =   -    (cos<t>    x      +    sintf>    x    )  sine       .  (4) 

Also  appearing  in  Eq.(l)  is 

r  /  2  ,_2    T.2 .  1/2      ,.     Tr2  „   2  .  2  .  ,_  . 

(w  /c   -  Kg)  7  ,     for   Kg  <  (j  /c  (5a) 

a0(u'¥  =        .  ,„2  2,    2,1/2     ,     „2 


2     2   2  1/2  2     2   2 

i(Kg  -  w  /c  )  7   ,   for   Kg  >  w  /c  (5b] 


The  amplitudes  of  the  incident  and  scattered  fields  are: 

*  <•.*«>  =  [*ll+*3  a  (A)  ^Bll  +  <*3  X  ^11  >  BX         <6> 

o v     II  ' 

and 

K"^ 
1S  (0,,%)  =  [%  -  x3  ^JTK^j-  ]An(W,%)  +  (x3  X  %)A1(o)/%)    (7) 

The  p-  and  s-  components  for  the  incident  wave  are  expressed  in 
terms  of   B|(   and   B.   ,  and  for  the  scattered  waves  in  terms  of 
the  Rayleigh  coefficients  A..   and  A. 

The  Rayleigh  coefficients  are  the  unknowns  in  the  problem. 
To  find  them  one  could  expand  the  field  in  the  metal  in  a  Rayleigh 
expansion  and  then  apply  the  boundary  conditions  at  the  grating 
surface  (this  is  the  Rayleigh  method) .   The  method  used  in  Ref .  8 
reduces  the  number  of  equations  in  half  by  employing  the  vectorial 
equivalent  of  the  Kirchhoff  integral  to  eliminate  the  field  in  the 
metal.   The  effective  boundary  condition  that  results  then  yields 
two  infinite  sets  of  linear  algebraic  equations  for  the  unknown 


coefficients.   These  are  Eqs. (2.18a)  and  (2.18b)  in  Ref.  8  and 
equations  (8a)  and  (8b)  in  Ref  9. 

The  equations  for  the  Rayleigh  coefficients  contain  the 
integral: 


I(«|G)  =  — -  SS    d  x„  exp[-iG-x..-iaC  (x  )  ] 


(8) 


c  cell 


where  a   is  the  area  of  a  unit  cell  on  the  bigrating  surface,  over 

which  the  integral  is  performed. 

The  perturbation  theory  begins  by  expanding  the  integrand  of 

"I"  to  first  order  in  C,  thus  yielding: 

I(a|G)  =  S^0  -  iaC(G)  (9) 

where 

C(G)  =  ~-  SS   d2x||C(x||)exp(-i^-x||)    .        (10) 
c  cell 

Under  the  assumption  that  the  dimensionless  parameter  af  is  small, 
we  can  substitute  equation  (9)  into  the  equations  for  the  A's,  and 
(assuming  that  f(0)=0)  find  the  first-order  results: 


A,i(3)  =is=-I'  n3"3')   [C32'V3'>  -a3*'Ai<3'>] 


a 


G' 


GG 


oa-± 

GG 


8a.o-lpa  c'(3) 

"3 


te2  Bira3,oBi] 


(lla) 


and 


a  (3)  =  loftY  c(3-3->   [ag3.A||(3')  +  bg3.A1(3-)] 

G' 


am 


g3>0-ip3^3) 
"3 


a3,o  Bn+b3,oBJ 


(lib) 


Here  the  prime  on  the  summation  indicates  that  G'?G,    and  we  have 
employed  the  notation  A..  (G)  for  A(|(&),K^)  (similarly  for  A.)  and 


1 


have  used  the  following  definitions: 

j3^  =  a(<o,Kg)  +  aQ((j/k||)  (12a) 

aGG'     =   a(w'K&     -a0(u'K2')  (12b) 

aG^  =  *3*(%  X  %'>  ^12c) 

bG^'  =  %  *  %•  (12d) 

KG  KG' 
CGG'  "  bG^'   +   a((o,K^)ao(W,K^.)  (12e) 

K-*  k 


a(o,Kg) 


2     .il/2 
e(co)-2  -  Kg 
c 


(12g) 


As  in  Ref.9,  we  assume  that  the  amplitudes  of  the  two 

components  of  the  specular  beam,  A..  (0)  and  A.  (0)  ,  can  dominate  the 

diffracted  orders.   Unlike  the  case  considered  in  Ref.9,  we  now 

consider  the  possibility  of  having  four,  not  two,  resonant 

evanescent  waves.   We  may  have  four  wavevectors  in  the  surface 

plane,  corresponding  to  four  different  reciprocal  lattice  vectors, 

which  satisfy  the  resonance  condition: 

%   =  ^11  +  Gr  =  f      (w)   ,        forr=l,...,4        (13) 
r  y 

where   K   (cj)   is  a  surface-polariton  wavevector  corresponding  to 

the  freguency  of  the  incident  light.   To  zeroth-order  in  the 

roughness  (the  flat  surface  limit),  the  locus  of  points   K   (w), 

sp 

for  a  given  freguency,  defines  a  circle.   Hence,  to  zeroth-order 
the  four  vectors   K-3   lie  on  a  circle,  at  the  four-fold  resonance. 

G 

If  any  one,  or  all  four,  of  the   K-3   are  at  resonance,  then  for 
that,  or  those,  wavevectors  we  have  to  zeroth-order: 


a(u,Kg  )  +  £(w)ao((J,Kg  )  =  0  (14) 

r  r 

We  are  therefore  going  to  assume  that  six  amplitudes  can  be 
larger  than  all  the  rest,  these  being   A..  (0)  ,  A.  (0),  A..  (G  )  , 
A..  (G  )  ,  A..  (G  )  ,  and  A..  (G  )  ;  and  for  each  of  these  we  can  write  an 
explicit  eguation,  using  Eq.(ll),  with  the  six  important  terms 
separated  out  from  the  sum  on  the  right  hand  side.   For  A.,  (1)  this 
procedure  gives  us: 


AU^    =  -^   i    [     I      f(l-p)clpA||(p)    +  r(l)[c10A„C0)    -      a1QA1(0)] 

p=2,3,4 
+   C(D[e1B||    -    a1QBi]    +   {    Y    ?  (1-j  )  [c^A,,  ( j  )  -aljA±  ( j  )  ] 

-  )   C(l-P)    a     A    (p)    ) 

p  J    J  (15) 


p=2,3,4 

To  simplify  the  notation  we  replace  the  variable  or  subscript   G. 

everywhere  by  j  alone  (where  j=l,..,4  are  the  resonant  beams  and 

j=0  the  specular  beam) .   The  prime  on  the  summation  indicates  to 

sum  on  all  j^0,l,2,3,4   .   The  nonresonant  terms  have  now  been 

segregated  in  the  brackets  on  the  right  hand  side.   The  equation 

for  A..  (2)  is  obtained  from  Eq.(15)  by  simply  interchanging  the  1 

and  2;  and  similarly  for  A..  (3)  and  A..  (4). 

The  specular  amplitudes  are: 
4 
Al<0>  =  ia00  1   n-r)a0rA||(r)  -  ^-B± 

r=l 

4 

+  ia00  {F  n-J)[aojA||(j)+bojAl(j)]  +  I      C(-r)b0rAl(r)}    (16 
i  r=l 


4 

lcoo  Z 
r= 

=1 

coo  ^Z 

f 

eoaoo  B 
pocoo   " 


+  ^{I  n-3)^0jAll^»-aOjAl(3)]  "I   f(-r)a0rAl(r)}  .  (17 


r=l 


The  amplitude  of  a  given  nonresonant  wave  can  also  be 

expressed  with  Eq.(ll),  but  now  we  drop  all  the  other  nonresonant 

terms  from  the  sum  on  the  right  hand  side: 

4 
A±(j)     =    ia^j    I    C(j-p)ajpAn(p)     +    f  ( j  )  [a..  0A„  (0)  +b^  QA±  (  0)  ] 
p=l 


A„(j)    = 


+   C(j)  [aj0B||+bj0B1] 


la  .  .      r    „ 

S^M    2    ?  CJ-p)cjpA||(p)    +    f(j)[cj0A„  (0)^^(0)] 


c 

"     p=l 


+   f  (j)  [ejB||-aj0B1]     }  (19) 


We  now  follow  the  development  given  in  Ref . (9) .   We 

substitute  Eqs.(18)  and  (19)  for  the  nonresonant  amplitudes  back 

into  Eq.(15)  for  the  the  resonant  amplitude  A(l) ;  and  for  the 

denominator   c    we  use 

K2  a(l)aQ(l)+K21 

Cll  =  1    +      o(l)ao(l)  =     a(l)ao(l)      '  (20) 

where   K.   now  stands  for   K^ .   and   «(j)   for   a(w,K-+.).   The 

denominator   c    vanishes  to  zeroth-order  at  the   K  -resonance. 

We  multiply  both  sides  of  the  resulting  equation  for  A..  (1)  by 

2     2      4 
a     (l)a     (l)-K-,   in  the  form 

a2(l)a^(l)-K^  =  [a(l)aQ(l)-K2]  [a ( 1) aQ ( 1 ) +K2 ]  (21) 

on  the  right  hand  side  and  as 
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2    2 

a2(l)aZQ(l)-Kl    =    e    ^  (  ~ 2 ~)        ■  (22) 

c     c 

on  the  left  hand  side.   Here  we  have 

wr  =  ^)c2Kr   '  (23) 

which  is  the  frequency  of  a  surf ace-polariton  at  wavevector   K 

on  the  flat  surface.   Finally,  we  note  that  in  the  equation  which 

we  now  have  for  A..  (1),   we  are  left  with  an  An(i)   on  the  right 

~2 

hand  side  in  a  term  proportional  to   f  ,   which  arose  from  the 

nonresonant  terms.   We  bring  this  term  over  to  the  left  side,  and 

use  it  to  renormalize  the  flat-surface  polariton  frequency: 

2  4 

~  2      2  f,   C  arr\      [  V   ?/    .  ? ,    ,      2 


03      =  CJ     <  1- 


r       (l+e)U2K  2  PP  rP 


I      r(r-p)C(p-r)appa] 


r     p=l 
(P^r) 


+  Y    ? (r-j)f (j-r)a 


:: 


c  .c  .     _ 
rj  jr.   2 

— J — ^— +a  . 

c  .  .    ri 

13 


for  r=l,2,3,  or  4.  (24) 

Here  we  employ  the  definition: 

Tr  =  a(r)aQ(r)  [a(r)oQ(r)  -   K^    ]  (25) 

Next  we  go  back  to  Eqs.(18)  and  (19)  for  the  nonresonant 
terms  and  substitute  them  into  Eqs.(16)  and  (17)  for   A. (0)   and 
A„(0)  . 

At  this  point  we  have  six  linear  algebraic  equations  for 
A..  (0)  ,  A.  (0)  ,  and  A.,  (r)  [r=l,  ...  ,4]   in  terms  of  B..  and  B .  ; 
expressed  in  matrix  form  these  equations  are: 
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2-2 

{a  -a     ) 

L12 

L13 

L14 

Ml 

Nl 

AH     U) 

L21 

2-2 

(«  -n2  ) 

L23 

L24 

M2 

N2 

A„(2) 

L31 

L3  2 

2-2 

(n  -o3z) 

L3  4 

M3 

N3 

A„(3) 

L41 

L4  2 

L4  3 

2-2 

(n  -n.   ) 

\             4     ' 

M4 

N4 

• 

AN(4) 

°1 

°2 

°3 

°4 

(1  +  P) 

Q 

A„(0) 

Rl 

R2 

R3 

R4 

S 

(1+T) 

A1(0) 

(Un-M,)B„  -  N^ 

"  N2B1 

"  N3B1 

-  N4B± 


v"l 

"1' 

<U2" 

-M2) 

Bll 

(u3- 

-M3) 

Bll 

<v 

■M4) 

Bll 

(W+V-P) 

Bll 

QB 


a 


00 


(X-S)B  -  (T+-^)B 
P0 


(26) 


In  this  equation  fi  =  w/(27rc/a)  =  a/X ;  Q  =  w  /(27rc/a)  ;  while  the 
other  symbols,  namely  the  capital  letters  L  through  X,  are  defined 
in  the  appendix. 

The  problem  now  is  to  solve,  numerically,  this  set  of  six 
equations  for  the  six  unknown  coefficients,  the  A's.   Then  these 
six  coefficients  can  be  substituted  back  into  Eqs.(18)  and  (19)  to 
find  the  nonresonant  coefficients.   With  all  the  coefficients  thus 
at  hand,  we  can  then  determine  the  scattered  field  vector- 
amplitude,   E  (K^) ,   from  Eq. (7)  —  for  each  diffracted  order  and 
for  each  evanescent  wave.   We  can  thus  obtain  the  reflectance,  or 
efficiency,  associated  with  each  of  the  grating's  diffracted 


12 


orders,  as  well  as  the  total  field  and  field-enhancement  near  the 
grating  surface. 

III.  Conclusion 

The  analytic  formulation  of  a  perturbation  theory  has  been 
presented  here,  in  complete  analogy  to  that  developed  in  Ref.9  , 
in  order  to  treat  the  special  case,  of  normal  incidence  on  the 
bigrating  surface,  which  was  not  treated  in  Ref.9  .   Numerical 
implementation  of  the  theory,  in  direct  analogy  to  that  in  Ref.9 
is  now  being  carried  out.   This  implementation  will  include  a 
comparison  of  the  perturbation  theory  to  the  exact  ,i.e., 
non-perturbative,  theory,  in  order  to  test  the  limits  of  validity 
of  the  former.   When  the  conditions  at  normal  incidence  are  such 
as  to  excite  simultaneously  four  surface-polaritons,  the 
absorptance  is  expected  to  be  enhanced  approximately  four-fold 
over  that  found  when  a  single  surface-polariton  is  resonantly 
excited. 

Again,  in  analogy  to  the  discussion  in  Ref.9,  the  results 
presented  here  can  be  used  to  calculate  the  dispersion  relation 
for  surface-polaritons  on  bigratings.   We  take  the  matrix  on  the 
left  hand  side  of  Eq. (26) ,  cross  out  the  fifth  and  sixth  rows  and 
columns  (corresponding  to  the  specular  beams) ,  take  the 
determinant  of  the  resulting  4x4  matrix,  and  set  it  equal  to  zero. 
The  primed  summation  symbols  that  appear  in  the  definitions  of  the 
matrix  elements  must  first  be  redefined  to  reinclude  the  specular 
beams,  j=0.   Solution  of  the  determinantal  equation  by  fixing   k.. 
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and  solving  for  complex  gj   will  then  yield  both  the  dispersion 
and  the  life-time  of  the  surf ace-polariton.   The  present 
formulation,  unlike  that  in  Ref.9,  will  enable  us  to  treat  those 
points  in  the   k..  -plane  where  two  Brillouin  zone  boundary  lines 
intersect,  that  is,  where  there  are  four  coupled  surface 
polaritons.   Work  on  numerically  solving  the  dispersion  relation 
is  just  getting  started. 
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Appendix 

Eq. (26)  employs  the  following  definitions  for  the  matrix 
elements. 


L        =  h 
rp  r 


if(r-p)c        -    )      C  (r-j )  C  (j-p)  (f    .c.    -a    .a.    )a. 

51  v      *'    rp        L      *  v      J '  *  KJ    *'  K    rj    jp      rj    jp'    jj_ 


(A-l) 


M 


=   h      if (r)c    .    +    )      f(r-p)f(p)a      a      a    _ 
r[    *  v    '    rO        L      *  v      *'*  Ktfl    pp   rp   pO 


p=l 
{*    r) 


~y      f(r-j)C(j)(f    -C    -a    .a..)a.. 


(A-2) 
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N      =    -h 
r  r 


iC(r)a    n    -   ]      C(r-p)f  (p)a      a      b    _ 


P=l 

(*    r) 


~y      C(r-j)f(j)(f    -a.n+a    .b.n)a..|  (A-3) 

A      * v      J/6  VJ/^    r]    ]0      rj    jO'    ]]J 


-a 


O      = 


00 


r  c 


00 


iC(-r)cQr   +   Y      C(P-r)C(-p)a      a 


PP   pr    Op 


p=l 

(*    r) 


~y      C  (j-r)f  (-j)  (fn  -c.    -a    .a. „)<*.. 


P   =   Sf    5    ^P)^-P)appaop   +   F    C(j)f(-j)(f0ic iQ+a    2)a 


(A-4) 


00 


a 


p=l 

4 


03    30       03'     33 


(A-5) 


00 


Rr   =    "a00 


—\    y    C(P)C(~P)a      a      b    _    +   y      C(J)C  (-j)  (frt.-b.„)a.rta.  . 

c       [   Z,    *  ^'^  v    ^'    pp   op   pO        L      *■  VJ/i>  v    J  '  v    Oj       jO'     30    33_ 


p=l 


(A-6) 


iC(-r)a.      -  y      C(p-r)C(-P)«      a     b   « 
j>  \      /    Or       L     *  V1^      /J>  v    ^'    pp  pr  pO 


p=l 

(*    r) 


~y      C (j-r)f(-j) (f •    an.+a.    b.rt) 
Z.         VJ       /iV,J/vjrOj      ]r    ]0 


a 


DD  J 


(A-7) 


S    =   a.J    y    C(P)C("P)a      a    nb    n    +   y      f  ( j  )  C  ("j  )  ( f  •  «-b  .  _)  a0  .a  .  . 
00|_    L    *>  V^/J>  v    ^'    pp   pO    pO         Z.      *"  ^J/i>  K    J'  v    jO      30'     O3    jj 

P=l  j 


T  =  °oo[  \  ^P)n-P)aPPbpo  +  F   f  «)f  M)  <  H^  +bjo  >ajjj 
P=i  j  d:i 

Ur  =  hr[iC(r)gr  -  J'    C(r-j)f (j)TjCrj 

j 
aoo  V'    ?,.,?,    -s 


V 


(A-8) 

(A-9) 
(A-10) 

(A-ll) 


00 
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w  = 


aooeo 

Coo^0 


X    = 


=  aocJ     ftJ>fM>Tjaoj 


(A-12 
(A-13 


In  the  above  equations  the  following  additional  definitions  have 
been  introduced: 

(A-14a) 


f   .  = 

raj 

c    . 

-      m3 

c  .  . 

33 

f .       = 

c  . 
-      lm 
c  .  . 
13 

gj  = 

cjo  - 

•   e . 
3 

h      = 
r 

~t    a 
r 

4 
c 
rr 

2 

ew 


27rc/a_ 


t  .  =  g.a  .  ./c  . 

D     :  DD  33 


(A-14b) 
(A-15) 

(A-16) 
(A-17) 


The  terms   a...  b..,  c..,  e.   were  defined  in  the  text,  Eq.(l2). 

They  are  dimensionless,  as  are   f . .  and  g.  .   The  summations  1' 
1  '  ID      ^3 

are,  as  in  the  text,  over  all  j  except  j=0, 1 , 2 , 3 , and  4  . 
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